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Abstract 

In contradistinction to a widespread belief, Wigner’s theorem allows the construction 
of true joint probabilities in phase space for distributions describing the object system as 
well as for distributions depending on the measurement apparatus. The fundamental role 
of Heisenberg’s uncertainty relations in Schrodinger form (including correlations) is pointed 
out for these two possible interpretations of joint probability distributions. E.g., in order 
that a multivariate normal probability distribution in phase space may correspond to a 
Wigner distribution of a pure or a mixed state, it is necessary and sufficient that Heisenberg’s 
uncertainty relation in Schrodinger form should be satisfied. 


1 Introduction 

Joint measurements of conjugate variables q and p are realized in many optical devices. This 
implies that one can think in this domain of a representation of quantum mechanics by means of 
joint probability distributions (j.p.d.) in the phase space of conjugate variables q and p [1]. This 
is perhaps the most convenient way to a realistic underpinning of quantum mechanics. A major 
advantage is that the incompatible variables q and p are c-numbers. The Wigner distribution 
function, which is widely used in optics, is the simplest language for coherent and squeezed states 
[2]. For these states the Wigner function is nonnegative. However, it is well known that the Wigner 
distribution cannot be considered as a true (nonnegative) probability distribution in general [3]. 
The aim of this paper is twofold: in the first part (sections 2 and 3) we present an analysis 
of the central question to consider phase space representations of quantum mechanics as true 
(nonnegative) probability distributions [4, 5] ; in the second part (sections 4 and 5) we emphasize 
the fundamental role of Heisenberg' & uncertainty relations in Schrodinger foim for Gaussian Wigner 
distributions and compare this with j.p.d. depending on the measurement arrangement (positive 
operator valued measures). 

2 Wigner’s theorem 

On account of the commutation relations between the operators q and p, there is no unique 
operator corresponding to the monomial q n p m . As a consequence there is no unique construction 
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of the j.p.d.. In general, a j.p.d. is completely determined by a given correspondence rule. 
Notwithstanding this arbit rareness, the existence of true probabilities in phase space is severely 
restricted by Winner’s theorem [3], which considers the following five requirements: 

(1) The j.p.d. is the mean value of an hermitian operator K(q } p ) depending on the c-numbers 
q and p: f(q,p) = tT[K(q } p)p]. 

(2) The j.p.d. is a linear functional of the density matrix (sesquilinear in the wavefunction): 
this means that K(q,p) is independent of p . 

(3) The j.p.d. is a true probability function: f(q } p) > 0. 

(4) When integrating over momentum space, the marginal distributions coincide with the 
proper quantummechanical probabilities in q: J f{q,p)dp =< q \ p \ q >. 

(5) When integrating over position, the marginal distributions coincide with the proper quan- 
tummechanical probabilities in p: / f{q,p)dq =< p | p | p >. 

Theorem 1 The five requirements (l)-(5) are incompatible . 

The requirement (2) is not explicitely present in the original version of Wigner’s theorem; the 
necessity of this requirement was emphasized by Mugur-Schlachter [6], who observed that in the 
absence of the arbitrary restriction (2) Wigner’s theorem cannot be realized. In the stronger 
version of Kruszynski and de Muynck [7] the requirement on one marginal distribution suffices. 

3 Realisation of positive phase space distributions 

For our purpose, it is sufficient to consider two different interpretations of j.p.d. as functionals of 
the density matrix. 

(1) The j.p.d. f(q,p) is interpreted as the probability that the variables q and p have certain 
values, the variable considered as a property possessed by the object system. In this case, two 
possibilities are left open for the construction of true j.p.d.: 

(1*1) w linear functional of p. 

In this case the requirements (l)-(5) are only compatible with a restricted class of functions. 
E.g. for the Weyl correspondence rule, the restricted class of functions are Gaussons (see section 4). 
The Wigner distribution cannot be considered as a true probability distribution in general, because 
e.g. it takes necessarily negative values for pure states that are not Gaussons. However, one 
can easily construct positive non-Gaussian Wigner j.p.d. corresponding to mixed states. For 
a representation of quantum mechanics by means of true Wigner j.p.d. one can add the new 
requirement that only nonnegative j.p.d. are physical states , This means e.g. that a one photon 
state is represented by a mixed state [10]. This idea is made plausible by the experimental fact 
that it is impossible to prepare a pure state with 100 % efficiency. 

(1.2) /(g,p) is a nonlinear functional of p. 

J.p.d. which are a nonlinear functional of the density matrix are not restricted by Wigner’s 
theorem. The j.p.d. which is the product of the proper quantum mechanical marginal distributions 
is a trivial example: f(q,p) =< q \ p \ q >< p \ p \ p >. Non-trivial examples with correlations 
exist also in the literature [11]. In this case the j.p.d. is a multilinear functional of the density 
matrix. We have considered a complete analysis of true distributions which are quadratic functional 
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of the density matrix [5]. This results in a new concept of j.p.d. which is based on a consistent 
phase space interpretation of the energy eigenstates of the wave function. 

(2) The j.p.d. /(?,p) is not function of the object system alone, but may also depend on the 
measurement arrangement of two incompatible observables Q and P . The measurements mutually 
influence each other in such a way that the singly measured quantum probability functions cannot 
be reproduced from the measurement results. In this case it is no longer desirable that the 
marginal probability distributions equal the single measured ones, hence Wigner’s theorem does 
not restrict this class of j.p.d. and f(q,p) may be a linear function of the density matrix. The 
optimal stochastic phase-space representations introduced by Prugovecki [12] are an example of 
this class. In general the distributions of class (2) can be considered in the framework of positive 
operator valued measures [13]. 

4 Heisenberg’s uncertainty relation in Schrodinger form 
and coherent and squeezed Wigner distributions 

We consider case (1.1) for the Weyl correspondence rule. In this case the construction of true 
j.p.d. for pure states is restricted by the remarkable and important theorem which was proven by 
Hudson [8] for one-dimensional systems and generalized by Soto and Claverie [9] for systems with 
an arbitrary number of degrees of freedom. 

Theorem 2 The necessary and sufficient condition for the Wigner distribution function of a pure 
state to be nonnegative is that the corresponding wave function < q\if> > is the exponential of a 
quadratic form . 

As a consequence the wave function represents a coherent or a squeezed state and the j.p.d. is 
a bivariate or a multivariate normal (Gaussian) distribution in phase space. Conversely, in two- 
dimensional phase space of the conjugate random variables q and p the most general normalised 
bivariate normal probability distribution with mean values q and p can be put in the standard 
form 

/( ) = ^^ exp { - 2Z ^ 9 ” ~ q K p ~ p ) + ~ P) 2 ] } 1 W 

where a q and a v and o qjt represent respectively the variances and the covariance a q = E\{q — g) 2 ] , 
etc.; E denotes the expectation value and A is the determinant of the covariance matrix: A = 
a q a p — cr^ > 0 . Schrodinger derived a more general and stronger form of Heisenberg’s uncertainty 
relation including the correlation a gj> : 


> ^ a /4 , (2) 

which we cadi “Heisenberg’s uncertainty relation in Schrodinger form”. It is easy to derive and 
to diagonalise the corresponding density matrix, f(q,p) may now represent a pure or a mixed 
state. The eigenfunctions < q\ip > are oscillator eigenfunctions functions multiplied by a common 
g-dependent phase factor which is characteristic for the correlation. We can show explicitely that 
there is a close connection between a Gaussian distribution in phase space, quadratic Hamiltonians 
and temperature dependent oscillator states. This implies a connection between physical and 
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statistical parameters. The eigenvalues of the corresponding density matrix are (1 - z)z n with 
z - (A - ft/2)/(A + k/2) , which leads to a sufficient condition for a bivariate normal probability 
distribution to be a quantum state: 


Theorem 3 In order that a bivariate normal probability distribution in phase space with variances 
o qi o p and covariance o qj> may correspond to a Wigner distribution of a pure or a mixed state, it is 
necessary and sufficient that Heisenberg’s uncertainty relation in Schrodinger form o q o p - o qiP 2 > 
h 2 /4 should be satisfied [4, 14 ] • 

It is very remarkable that the Schrodinger form of Heisenberg’s uncertainty relation, which is a 
necessary condition to be fulfilled for every Wigner distribution function, is also a sufficient condi- 
tion in the case of a bivariate normal probability distribution. Indeed, to be a Wigner distribution 
function, f(q,p ) must satisfy an infinite set of KLM [15] or equivalent conditions in general, but 
for the two-dimensional Gaussian distribution the infinite set reduces to one simple necessary and 
sufficient physical condition. In this respect, the uncertainty relation in Schrodinger form is more 
fundamental than Heisenberg’s relation in the usual, less stronger form o q cr p > h /4. Moreover, 
the Schrodinger form is invariant for linear canonical transformations (in general Sp{2n, R) in- 
variant transformations), while the usual form is not. Finally, for quadratic Hamiltonians, which 
are closely related to the Gaussian Wigner distribution, the Schrodinger form remains invariant 
during the motion if the variances and the covariance are dependent on time. Indeed, in this case 
the quantum Liouville equation is equivalent to the classical Liouville equation and therefore q, 
p,o q> Op and &qj> have the same time dependence as in the classical case. These are further reasons 
why the uncertainty relation in Schrodinger form is more relevant than Heisenberg s relation in 
the usual form. 

For systems with an arbitrary number of degrees of freedom the strong form of Heisenberg s 
uncertainty relation is derived from the inequality tr(o f pa) > 0 where the vector a is given by 
a = A(q — q) + B(p — p), A and B being arbitrary matrices, and which takes the form: 




o q>q o qj) ih/2 

Op, q — lh/2 Opj, 


A 

B 


> 0 . 


(3) 


Therefore Heisenberg’s uncertainty relation in Schrodinger form takes now the matrix form: 


o - thP/2 > 0. 


(4) 


where o is the covariance matrix I '*■« I and ft the fundamental symplectic matrix I ‘I- 

1 9 P V P*P I I ' 

Theorem 4 The necessary and sufficient conditions for a Gaussian phase space function to be a 
Wigner distribution is that the covariance matrix o satisfies Heisenberg’s uncertainty relation in 
Schrodinger form: o - ihfi/2 > 0 [4]- 

Analogous remarks as for the bivariate j.p.d. are valid for the multivariate j.p.d., the eq. 4 is 
now Sp(2n, R) invariant. The theorem entails a considerable simplification with respect to the 
theorem of Simon, Sudarshan and Mukanda [17], where Sp(2n, R) invariant powers of (3cr satisfy 
n complicated inequalities. The difference between a pure and a mixed state is given by a theorem 
of Littlejohn [16]: 

Theorem 5 The necessary and sufficient condition for a Gaussian Wigner distribution to be a 
pure state is that the matrix 2 o/h is a symplectic matrix: oflo = (h 2 /4){3. 

In two dimensions the matrix relation reduces to o q o p — o q # = h 2 / 4. 
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5 Heisenberg’s uncertainty relation in Schrodinger form 
for j.p.d. depending on the measurement arrangement 

It was argued in section 3 that the construction of j.p.d. of class (2) is not restricted by Wigner’s 
theorem. Requiring Galilei invariance, linearity and positivity for any density matrix describing 
the object system, we have for the most general form of the j.p.d.: 

/(?> P) = ^ tl(Dq J> p rneat Dqj,p o l >] ) (5) 

where p mea » and pdj are the density matrices describing exhaustively the measurement apparatus 
and the object system and D qj) represents the displacement operator. If both Pmea$ and p&j are 
pure states then f(q,p ) reduces to the transition probability f(q,p) = ft _n tr( D qj> yp ,*,])■ The 
marginal distributions are always given by the convolution of two true probability densities: 

J /(g»p)dp =< q I Pm**, I q > * < q | pcbj | q >, (6) 

J /(g.p)dg =< p I Pmea, \ P > * < P | pcbj \ P >, (7) 

which can be seen as accuracy calibrations given by the confidence functions < q \ pmeaM | 9 > 
and < p | pmea# | P >• The couple q,< q | Pmea# | 9 > together with p, < p | Pmea# I P > can also be 
interpreted as a fuzzy sample point in phase space [12]. Remark also that, for these j.p.d. the 
ordening of operators is equivalent with a measuring procedure. One can also write the j.p.d. as 
a convolution of two Wigner distributions: 

/(9>P) = frnea§{<l}P) * fcbj(9)P') 1 (®) 

the first one representing the measurement procedure and the second one describing the object sys- 
tem. This “smoothing” or “coarse graining” of the Wigner distribution eleminates fast oscillations 
in h and gives therefore a better representation in the classical limit [18]. Another consequence 
of the last formula is that the covariance matrix a is the sum of the covariance matrix of the 
object system and the of the measurement procedure. Hence we obtain the “operational” 

uncertainty relation 

a — ih/3 > 0. (9) 

which reduces in one dimension to a q a p — a qj> 2 > h 2 . This operational uncertainty relation is 
in accordance with the experimental uncertainty relation (Ag)ex(Ap) cx ~ h [19]. Comparing 
this with the uncertainty relations for the j.p.d. of the preceding section, we observe that the 
inequalities are the same, except for the essential difference that ft replaced by 2ft, expressing the 
presence of extra noise due to the measurement procedure. 
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